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Abstract. 

We study the phase-ordering kinetics following a temperature quench of 0{M) 
continuous symmetry models with A/" = 3 and A/" = 4 on graphs. By means of extensive 
simulations, we show that the global pattern of scaling behaviors is analogous to the 
one found on usual lattices. The exponent for the integrated response function and 
the exponent z, describing the growing length, are related to the large scale topology 
of the networks through the spectral dimension and the fractal dimension alone, by 
means of the same expressions provided by the analytic solution of the J\f ^ oo limit. 
This suggests that the large- A^ value of these exponents could be exact for every TV > 2. 
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1. Introduction 

A ferromagnetic system quenched from a high temperature disordered state to an 
ordered phase with broken ergodicity evolves via phase-ordering dynamics. In the late 
stage of growth, the correlation of the order parameter between sites i,j at times s,t 
can be expressed as the sum of two terms P 

C.,{t,s) = Ct^{t-s)+Ctf{t,s). (1) 

The first term describes the contribution provided by degrees of freedom which in 
the interval [s, t] are not interested by non-equilibrium effects, such as the passage 
of topological defects, while the second contains the non-equilibrium information. 
Analogously, also the integrated response function, or zero field cooled magnetization, 
measured on site i at time t after a perturbation has been switched on in j from time s 
onwards, takes an analogous addictive form [2] 

X,,it,s)=x^it-s) + x-!it,s). (2) 

On regular lattices, due to space homogeneity and isotropy, correlation and response 
function depend only on the distance r between i and j. One has, therefore, Cij{t, s) = 
C{r, t, s), and similarly for Xijity -s)- 

The non-equilibrium behaviour is characterized by the dynamical scaling symmetry, 
a self-similarity where time acts as a length rescaling. When scaling holds, the states 
sequentially probed by the system are statistically equivalent provided lengths are 
measured in units of the characteristic length scale L{t), which increases in time. All the 
time dependence enters through L(t), and the aging parts in Eqs. ( |T|2|) take a scaling 
form in terms of rescaled variables llj x = r/L{s) and y = L{t)/L{s) 



r^{r,t,s) = C{x,y), (3) 

X^'ir,t,s) = s-''-xix,y). (4) 

The characteristic length usually grows according to a power law 

L(t)~tV^ (5) 



Interestingly, non-equilibrium exponents are expected to be universal, namely to depend 
only on a restricted set of parameters. On regular lattices, where a substantial 
understanding of the dynamics has been achieved by means of exact solutions, 
approximate theories and numerical simulations [3], these exponents depend only on 
the space dimensionality, the number of components JV of the order parameter and the 
conservation laws of the dynamics. In this paper we will consider primarily z and a^. It 
is known [3] , that z depends only on the conservation laws, being z = 2 in the case of a 
non conserved order parameter considered in this paper. Regarding a^, for continuous 
symmetry models (A/" > 1) with non-conserved order parameter, it was conjectured in 
|4J to obey 

for d<A 

1 with log corrections for d = A (6) 
1 for d> 4. 
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Let us notice that these exponents share the property of being TV-independent. Their 
value therefore can be computed in the soluble large- A/" model |5]. The reference 
framework of a soluble theory is a great advantage in the study of phase-ordering on 
networks or inhomogeneous graphs considered in this paper. In fact our understanding 
of these systems is largely incomplete, although examples can be found in disordered 
materials, percolation clusters, glasses, polymers, and bio-molecules, and are also present 
in interdisciplinary studies, ranging from economics to chemistry and social sciences [7j. 
The aim of this paper is to investigate if also on networks the non-equilibrium dynamics 
of statistical models takes a scaling structure, and its universality features. Moreover, 
it is interesting to understand which topological indices of the graph play the role of the 
euclidean dimension d on regular lattices in determining universal quantities. We will 
restricted our attention to physical graphs [8] : These are networks with the appropriate 
topological features to represent real physical structures, namely they are embeddable 
in a finite dimensional space and have bounded degree. 

Equilibrium properties of models defined on physical graphs, and in particular 
the relevance of their topology, are quite well understood. As far as systems with 
continuous symmetry are concerned, such as 0{Af) models {Af > 2), a unique parameter, 
the "spectral dimension" ds, encodes the relevant large scale topological features of 
the network and regulates the critical properties, dg is related to the low eigenvalues 
behaviour of the density of state of the Laplacian operator [9] , and can be considered for 
these models as the topological indicator replacing the Euclidean dimension d on graphs: 
The spectral dimension univocally determines the existence of phase transitions |10j and 
controls critical behaviours [llj, much in the same way as the Euclidean dimension d 
does on usual translation invariant lattices. 

Regarding non-equilibrium, in the large- A/" model it was shown fl2[ [13] that the 
general framework of scaling behaviour discussed above is maintained on generic graphs. 
In particular the topology of the network enters the exponent z only through the fractal 
dimension df and the spectral dimension ds, 

z = 2df/ds, (7) 

while depends only on the spectral dimension obeying ([6]) with ds occurring in place 
of d. 

In this paper we want to complement the large- A/" analysis by studying the phase- 
ordering kinetics of some 0{N') vector models with finite AT on a class of physical graphs. 
We will consider geometrical fractals without phase transition at finite temperature, 
such as the Sierpinski gasket or the T-fractal, and others structures, obtained by direct 
products among graphs |8j, which on the contrary feature a phase transitions at a finite 
temperature T^. Our simulations of these systems evolving with relaxational dynamics 
(non-conserved order parameter) show that, for A/" = 3 and A/" = 4, the dynamical 
exponents z and are correctly predicted by the large- A/" model and then depend only 
on the fractal and spectral dimensions df,ds. This suggests that in presence of scaling 
these exponents take the same value for any A/" > 2, much in the same way as it happens 
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on euclidean lattices. 

This paper is organized as follows: In Sec. [2]we introduce the physical graphs giving 
a definition of fractal and spectral dimension. In Sec. [3] we define the 0{N') models 
that will be considered in the simulations. We also introduce the basic observables, 
and discuss the numerical techniques. In Sec. H] we present our results for different 
structures. Sec. [5] contains a final discussion and the conclusions. 

2. Physical Graphs 

A graph (network) ^ is a discrete structure defined by a set of sites i connected pairwise 
by unoriented links (edges) {«, j}. The chemical distance r^j-, i.e. the number of links 
in the shortest path connecting sites i and j, naturally defines on ^ a metric. Van 
Hove spheres , allowing to explore large scales of the graph, can be constructed using 
this metric. The van Hove sphere So^r of radius r and center a is the sub-graph of Q 
composed by the sites whose distance from o is smaller than r + 1. Calling No^r the 
number of sites in So^n the fractal dimension of the graph is defined by its asymptotic 
behaviour at large scales 

No,r ~ r'^f (8) 

where ~ denotes the behaviour for large r. In the following we will consider only physical 
graphs embeddable in a finite dimensional space {df is well defined and finite). Moreover 
we require that the degree Zi (number of neighbours of the site i) is bounded. 

Algebraic graph theory provides powerful tools for the description of the topology 
of generic networks, by means of characteristic matrices. The adjacency matrix Aij of 
a graph has entries equal to 1 if z and j are neighbouring sites ({«,j} is a link) and 
Ai^j = otherwise. The Laplacian matrix Aij is defined as 

where Zi = J2j is the degree of i. Interestingly, Aij is the generalization to graphs 
of the usual Laplacian operator of Euclidean structures [8j. In particular its spectrum 
is positive and the constant vector is the only eigenvector of eigenvalue zero. Moreover, 
on physical graphs the spectral density p(/) of Aj^-, is expected to behave as [9] 

p(/) ~ /'^^/^-i ^10) 

where ~ denotes the behaviour for small Vs. Eq. ( JTOl) defines the spectral dimension ds 
of the physical graph. 

In order to build graphs of larger dimensions we introduce the direct product 
between graphs. Given Q and H the direct product ^ x 7i is a graph whose sites 
are labelled by a pair (i, j) with i and j belonging to Q and H respectively. {i,j) and 
are neighbour sites in Q x Ti, if i = i' and {j,j'} is a link of 7i, or if j = j' and 
{i, i'} is a link of Q. A basic properties of ^ x H is that, calling d^ and d^ the dimensions 
of the graph Q, one has [8] 

dS-^ = d^,+d^ df'^ = d'^. + d'} (11) 
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T-FRACTAL SIERPINSKI-GASKET 



Figure 1. The finitely ramified fractals considered in the paper. 

i.e. the spectral and fractal dimensions of the product graph are the sum of the 
dimensions of the original graphs. 

In the simulations, we will consider models defined on graphs with known dg 
and df in order to verify the relevance of these dimensions, describing the large 
scale topology, on phase ordering. In particular, we focus on two finitely ramified 
fractals [H] with dg < 2, the T-fractal and the Sierpinski gasket, whose fractal and 
spectral dimensions can be analytically evaluated by means of exact renormalizations 
m, yielding df = log(3)/ log(2), 4 = log(9)/ log(6) and df = log(3)/ log(2), 4 = 
log(9)/log(5) respectively. Moreover, in order to explore the region between d = 2 and 
d = 3, we also consider the graphs obtained from the product of two gaskets and two 
T-fractals, which both have spectral dimensions 2 < ds < 3 . 

3. The Model 

The Hamiltonian of 0{N') spin systems on a graph is given by: 

N 

Ho[a] = -JY^ Ai^jdi ■ Of (12) 

where |crj| = 1 are unitary A/'-component vectors, and is the number of sites in the 
network. In the following we will set J = 1. Phase ordering is obtained by evolving 
an initially fully disordered configuration with a dynamics at a temperature T where 
the system presents an ordered phase at equilibrium. In the following we will consider 
quenches to T = 0. This guarantees that the dynamics is of the phase-ordering type even 
for those systems with (is < 2, namely without phase transition at finite temperature. 
Phase ordering can be characterized by the correlation function, defined as 

C,,(t, s) = {a,{t) ■ a,{s)) - {a,it)) ■ (13) 

Recalling Eqs. f|T][3|) . the non-equilibrium scaling properties are encoded in the aging 
part C°^(t, s) which, in turn, can be obtained by subtracting the equilibrium contribution 
C-j{t,s) from the whole correlation Cij(t,s). Notice however that at T = 0, C-j(t,s) 
vanishes and hence C^jit, s) = Cij{t, s). We will drop therefore the superscript ag from 
the correlation function in the following. In order to simplify the analysis one usually 
restricts the attention to Cij{t,t), namely the equal time correlation function, or to 
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the autocorrelation function Cii(t,s). On lattices Cij(t,t) = G{r,t) is a function of the 
euclidean distance r between i and j which, according to Eq. ([3]) scales as 

G{r,t)=g{r/L{t)) (14) 

This allows one to extract L{t) as, for instance, the half-height width of G{r,t). 
Concerning the autocorrelation function, on lattices it does not depend on i, due to 
translational invariance, so that Cii{t,s) = C{t,s) with, following Eq. ([3]), the scaling 
form 

C{t,s) = f{y)r^y-\ (15) 

where A depends both on d and on A/" [3] . 

When phase-ordering occurs on graphs one has the additional feature of the 
dependence of Cij{t,t) on i and j (and Cuit.s) on i), complicating the analysis and 
hindering the scaling properties. Therefore, in order to simplify the study, we will 
resort in the following to particular correlations, with a transparent physical meaning, 
which where shown in p!3t [T6] to be useful in detecting scaling properties. More 
precisely, for the equal time correlation function, on the Sierpinski gasket and on the 
T-fractal we will compute it restricting i and j on the baseline of the structure (see 
Fig. [1]), as this procedure allows to soften the log-periodic oscillations characteristic of 
deterministic fractals [I7j and gives the best results for the scaling plots, as will be shown 
in Sec. HI Regarding the autocorrelation function we will consider its spatial average 
Jj Sill ^m(^) ■s), in order to obtain the best statistic. For simplicity these quantities will 
be denoted as G{r,t) and G(t, s) in analogy to their counterparts on lattices. 

The response of the system to an external field can be studied introducing the 
susceptibility 

Xij{t, s) = \im ^{ai{t))h^(^,) ■ Hj (16) 

where {(Ti{t))h^(^s) denotes the expectation value at time t of the spin cxj when a field 
hj = hfij {fij being a unitary vector), changing the Hamiltonian to H[a] = HqIct] —aj ■ hj 
is switched on from time s < t onwards on site j. Differently from what happens for 
Cij(t,s), the equilibrium contribution to Xtjit^s) does not vanish. Then, in order to 
isolate the scaling part Xijity s), according to Eq. ([2]) one has to subtract Xij(t~s) from 
the whole response measured during the quench. x|j(t — s) can be numerically evaluated 
as the response of a system prepared initially in the equilibrium state at T = 0, namely 
with all the spins aligned. Actually, if one does this the system starts to evolve as soon 
as a small field /i > is switched on even if T = 0, originating a non-vanishing response. 
Considering, for simphcity, the autoresponse x(^5 s) = Xaity ■s), on lattices, from Eq. (jlj) 
one finds 

x'''{t,s) = s-''-h{y) (17) 

with depending only on d according to Eq. ([6]). On networks the response function 
is site dependent, so, similarly to what done for the autocorrelation function, we will 
consider the spatially averaged quantity jj^ X^ili Xaity ^) which will be denoted as x{t, s). 
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Zero temperature dynamics can be implemented in different ways. For example one 
can set T = in a Metropolis updating rule. A convenient choice adopted in this paper 
consists in aligning a randomly chosen spin cij with the local field, i.e. the spin ai is 
turned into 



With this choice, at each move the local energy is minimized. Such a simple updating 
rule proves to be very efficient, on networks and lattices as well. For instance it allows to 
reproduce easily the analytically known behaviour of the 0(2) model in one-dimension 
|18j . whereas the Metropolis rule fails due to extremely long transients (see Appendix 
for details). Notice that in one- dimensional 0(2) model scaling is violated, and our 
hypothesis on dynamical exponents does not hold. 

The scope of this study being the analysis of the behaviour of finite J\f models, 
the less numerically demanding case to start with would naturally be the 0(2) model. 
However, as discussed in the Appendix, the dynamics of this model is pinned on self- 
similar graphs such as the Sierpinski gasket or the T-fractal considered here. Then, 
in order to study the phase ordering kinetics in this case, one should allow the system 
to depin by quenching to a small but finite temperature. Since the above mentioned 
structures, however, are disordered at any finite temperature, in order to recover true 
phase-ordering one should then let T — > 0, in the same way as for phase-ordering in 
the one-dimensional Ising model with Kawasaki dynamics [19]. Since this procedure is 
very numerically demanding, hereafter we will focus on the next simplest cases, namely 
with A/" = 3 and A/" = 4 (the latter restricted to the Sierpinski gasket). Performing 
extensive numerical simulation of these models on different self-similar graphs, we show 
that phase ordering exhibits dynamical scaling, similarly to what is known on lattices, 
in any case. In particular, the correlation length grows according to Eq. ([5]) and the 
two time correlation and response functions scale as ( JT5i) and ( IT71) . We also find that 
the exponents z and are well consistent with their large- A/" value. This suggests that 
also on graphs, as on lattices, these exponent do not depend on A/", and that their value 
can be predicted by the solution of the large- A/" model for every scaling system with 
Ar> 2. 

4. Results 

In the numerical simulations we considered the Sierpinski gasket, the T-fractal, and the 
graphs obtained by a direct product of two gaskets and of two T-fractals. The number 
of sites of the structures is 2391486, 1549324,10771524 and 4787344 respectively. A 
time-step is made of A^ elementary Montecarlo moves ( JTSl) . The data are obtained by 
averaging over at least 500 dynamical realizations, each with a different random initial 
condition. For the calculus of xi^yt) ^ different random field Ui is applied in each run. 
To evaluate the limit /i — in Eq. (fT6l) . we have simulated systems with small external 




(18) 
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Figure 2. Scaling of G{r, t) evaluated on the baseline of a Sierpinski gasket. 



fields h and then we have verified that s) is independent of h (e.g. comparing 
systems with field h and h/2). In our simulations h varies between 10~^ and 10~^. 

In Figure [2] we show the behaviour of the equal time correlation function G{r, t) 
defined in Sec. [3] for the Sierpinski gasket. The Figure evidences an excellent scaling, in 
the form expected from Eq. f|T^ . Quite surprisingly, the scaling is very good even for 
very short times {t ranges from 8 to 46600) evidencing the efficiency of the dynamics 
considered. Analogous results were obtained for the T-fractal. In Figure [3] we plot the 
growth of L{t) as a function of time, finding a very good agreement with the asymptotic 
behaviour L{t) ~ fd.s/(2df) (^Yiamelj t^-^^i g^j^^ ^o.ss? ^^ie Sierpinski gasket and the T- 
fractal respectively). We note that superimposed to the expected behaviour there are 
small log-periodic oscillations, typical of fractal structures [I7j. If the average of the 
correlation function is not restricted on the baseline of the structure, but is taken on the 
whole fractal, the log-periodic oscillations are more evident, as shown in Figure[3]for the 
Sierpinski gasket, while keeping the same slope for the L[t). If this average procedure 
is used, it is more difficult to extract the scaling behaviour of the correlation function. 

Let us now turn to the response function. It has been analyzed as follows: After 
computing x"'^{'t:^) by subtracting s) to the full non-equilibrium response, as 

explained in Sec. [3l we fix the ratio y = t/s and plot x'^^i't, s) versus s (Fig. Hj). In doing 
so we obtain, from f|T7|) . an estimate for from the slope of the plot. As illustrated in 
Fig. m both for the Sierpinski gasket and for the T-fractal the value of is consistent 
with ds/2 — 1 (namely = —0.317 for the Sierpinski gasket and = —0.387 for the T- 
fractal) i.e. the analytic expression obtained in the A/" — >■ oo case (except for the smaller 
values of t/s, as we will discuss below). Finally, we verified the scaling relation (ITTl) by 
plotting ■s) versus t/s for different values of s and checking for data collapse, as 

shown for the Sierpinski gasket in Fig. [5l The collapse is good for t/s > 5, whereas it is 
poor for smaller values of t/s. This can be interpreted as due to finite-s preasymptotic 
corrections, which are more effective at small t/s, similarly to what observed on lattices 
[H E] . Analogous result can be obtained for the T-fractal. Let us stress the particular 
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Figure 3. Growth of L(t) as a function of time, anomalous diffusive behaviour 
L{t) oc t-^/^ is well verified. For distances evaluated along the baseline, best fits yield 
1/z ^ 0.42 ± 0.01 and 1/z = 0.39 ± 0.01 for the Sierpinski gasket and the T-fractal, 
respectively, very well consistent with the large-TV behaviour 1/z — dg/{2df) in both 
cases. For the Sierpinski Gasket we also plot L{t), evaluated averaging over sites 
belonging to the bulk of the structure. The value of the exponent does not change, 
but the log-periodic oscillations are much larger. 

feature of the cases with (i^ < 2 of x^^(t—s) (and hence ■§)) scahng itself as s~"'^/i(t/s), 
with the same exponent of x'^^i^-, and with a scahng function h{y) that, ahhough 
different from h{y), has the same asymptotic behaviour h{y) ~ y'*^^ for y ^ 1. In this 
case, therefore, one gets the same information on the exponents from s), x'^^i^iS) 
and ~ The situation is different for ds > 2. Here the stationary term converges 
to the finite equihbrium susceptibihty, hm^^^oo X^^i^ ~ ^) = Xeq, whereas s) scales 

with a positive (see below). In this case, therefore, xit^^) is dominated by x'^^i't/s) 
and the subtraction of x^*(V"^) is necessary in order to make the scaling of the aging 
part manifest. In Fig. [6] we evaluate the exponent a^, following the procedure discussed 
above, both for the product of gaskets and for the product of T-fractals. Also in this 
case we get a very good agreement with the expected behaviour = ds/2 — 1 (namely 
= 0.365 and = 0.226 for the product of gaskets and the product of T-fractals 
respectively) . 

In order to provide a further check on dynamical scaling, we have considered also 
the two-time correlation function C{t, s) introduced in Sec. [3l The scaling relation (ITSjl 
is well verified, as shown in Fig. [7|for the Sierpinski gasket. A residual correction, due 
to the finite values of t^ can be observed for large t/tw, similarly to what is known 
in lattices [20]. The quality of the scaling improves by increasing t^, with the curves 
for tui = 40, 80, 160 almost collapsing, as expected. Similar results were found for the 
T-fractal. 

Up to now, by focusing on the case N" = 3, we have shown that the global pattern 
of behaviours observed on lattices is maintained on the fractal structures we have 
considered, both with dg < 2 and with dg > 2. Moreover, exponents such as z and 
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Figure 4. Estimate of the exponent for the Sierpinski gasket and the T-fractal (see 
text). Best fits yield = —0.32 ±0.01 for the Sierpinski gasket and a-^ = —0.38 ±0.01 
for the T-fractal. In both cases these values are consistent with the hypothesis 
— dg /2 ~ 1. 




Figure 5. Scaling of the function x'^^i'tjs) on the Sierpinski gasket. At large time 
ratio, a good agreement with the scaling hypothesis pT]) is verified for s > 10. 



that are known on lattices not to depend on Af, are found consistent with the large- 
Af value. These results suggest that the property of being A/"-independent may hold for 
any scaling system with M >2. To check this hypothesis, we have considered the 0(4) 
model on the less computational demanding structure, i.e. the Sierpinski gasket. In Fig. 
E] we show the growth of the correlation length, which is again in good agreement with 
an exponent z = 2df/ds. Analogously, also x°^(^,s), at large enough time ratios, scales 
as ( flTl) with an exponent very well consistent with the large- A/" result = ds/2 — 1, a.s 
shown in Fig. [9l 
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Figure 6. Estimate of the exponent for the product of Sierpinski gaskets and the 
product of T-fractals. Best fits yield = 0.36 ± 0.01 for the product of Sierpinski 
gaskets and — 0.22 ± 0.01 for the product of T-fractals. In both cases a good 
agreement with the large-A/" value is evidenced. Notice that in the first case the 
asymptotic regime is reached much earlier than in the second one. 




Figure 7. Scaling of the two time correlation function C{s,t). A good agreement 
with the hypothesis C{s,t) = h{s/t) is found. 



5. Conclusions 

In this paper, we have studied the phase-ordering kinetics of 0(3) and 0(4) models on 
self-similar physical graphs with known spectral dimension. On homogeneous lattices 
these models obey dynamical scaling. We have shown by means of extensive numerical 
simulations that the same symmetry appears to be obeyed also on the networks 
considered. The non-equilibrium exponents z and a^^ are consistent with the value 
provided by the large-A/" model. These results suggest that Equations (El) and (171) for 
z and may hold for any Af > 2, for systems interested by dynamical scaling. It 
would important to substantiate this hypothesis by means of analytic calculations that 
go beyond the large-A/" limit by accessing directly the case of finite A/". This could be 
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Figure 8. Growth of L{t) as a function of time for the 0(4) model on the Sierpinski 
gasket, a comparison with the expected behaviour cf^^/C^'^f') is provided. 




Figure 9. Scaling of the susceptibility x°'^itis), for large t/s a good agreement with 
the scaling hypothesis PT|) with = (is/2 — 1 is evidenced. 



possibly achieved by using a Gaussian Auxiliary Field approximation on the equation 
of motion for the 0{Af) model [2T] . 

On the other hand, concerning the issue of the generality of the scaling property, 
it would be interesting to study if the breakdown of dynamical scaling always occurs in 
the 0(2) model on structures with dg < 2, similarly to what observed on homogeneous 
one-dimensional lattices [T8] . 

Appendix A. 0(2) model 

In this appendix we present some results obtained for the 0(2) model using zero 
temperature heath bath dynamics (fTSl) . First we discuss phase ordering in a one- 
dimensional system. We show that our dynamics reproduce the asymptotic behaviour 
predicted analytically in [18] where zero-temperature phase ordering for the 0(2) model 
is studied by means of continuous time dissipative dynamics. In [18] it is evidenced that 
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Figure Al. The growth of the correlation lengths L{t) and LK{t). 



the dynamical scaling is violated since two different correlation lengths are present. The 
first one, the phase winding length L{t), is the characteristic length of the spin spin-spin 
correlation function defined by Eq. (fT^ . The other is the phase coherence length Lx{t) 
representing the typical distance for which the winding of the spins along the system 
changes its direction. In particular, calling — vr < 6i < n the angle between neighbouring 
spins (jj and aj+i, one defines on each site the binary variable ki = ±1 if 6*4 > or 6'j < 
respectively. According to llSj, the correlation function GK{r,t) = {ki(t)ki+r(t)) scales 
differently from G{r,t); in particular 

G{r,t) = gir/m) = g{x) G,,(r,t) = gk{r/LK{t)) = gu{x) (A.l) 

with L{t) ~ and Lx(t) ~ t^/'^. In our simulation with zero temperature heath bath 
dynamics we verify the laws ( lA.ll) . In particular in Fig. lAll we display the growth in 
time of the two correlation function evidencing a very good agreement with the expected 
results. We remark that adopting different dynamics, e.g. Metropolis, the analytically 
known behaviour is not easily reproduced in simulations because of very-long transients 
(only the asymptotic behaviour is expected to be independent of the transition rates). 

As explained in Sec. [3] we have initially considered the 0(2) rotator model on 
fractals, and we obtain that the system gets frozen into metastable states. In particular 
the energy E{t) does not tend to zero and the correlation length does not diverge. Figure 
A2I shows this behaviour for the Sierpinski gasket. Analogous behaviours are present on 
other fractals. 
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